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We investigate the influence of a dissipative environment which effectively comprises the effects of 
counterions and hydration shells, on the transport properties of short DNA wires. Their electronic 
structure is captured by a tight-binding model which is embedded in a bath consisting of a collec- 
tion of harmonic oscillators. Without coupling to the bath a temperature independent gap opens in 
the electronic spectrum. Upon allowing for electron-bath interaction the gap becomes temperature 
dependent. It increases with temperature in the weak-coupling limit to the bath degrees of freedom. 
In the strong-coupling regime a bath-induced pseudo-gap is formed. As a result, a crossover from 
tunneling to activated behavior in the low-voltage region of the I-V characteristics is observed with 
increasing temperature. The temperature dependence of the transmission near the Fermi energy, 
t(EF), manifests an Arrhenius-like behavior in agreement with recent transport experiments. More- 
over, t(Ep) shows a weak exponential dependence on the wire length, typical of strong incoherent 
transport. Disorder effects smear the electronic bands, but do not appreciably affect the pseudo-gap 
formation. 
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I. INTRODUCTION 



The idea that conduction pathways in DNA molecules 
may be built up as a result of the hybridization of the 7r 
orbital stack along consecutive base pairs can be traced 
back to the 1960'silt was not, however, till recently that 
a revival of interest on DNA as a potential conductor oc- 
curred. This was mainly triggered by the observation of 
long-range electron transfer between intercalated donor 
and acceptor centers in DNA molecules in solution. 2 
Subsequent experimental results 3 ^^^*^^^ were contro- 
versial as they showed different functional dependences 
of electron transfer rates on the donor-acceptor separa- 
tion. Thus, strong exponential fall-off 4 ^ typical for su- 
perexchange mediated transfer as well as a weak, alge- 
braic dependence 3 ^ characteristic of sequential hopping 
processes were reported. Meanwhile, theoretical work 
has led to an emerging picture where different mecha- 
nisms may coexist depending on base-pair sequence and 
energetics liSiii 

In parallel to these developments in the chemical 
physics community, DC transport experiments on 
A-DNA as well as on poly(dG)-poly(dC) and poly(dA)- 
poly(dT) molecules between metal electrodes have been 
performed, 12 ! 13 : 14 : 15 : 16 : 17 : 18 : 19 Several fundamental diffi- 
culties have to be surmounted in this kind of experi- 
ments: (i) how to create good contacts to the metal 
electrodes, (ii) how to control charge injection into the 
molecule, (iii) single molecule vs. bundles of molecules 
and (iv) dry vs. aqueous environments, among others. 
Consequently, sample preparation and the specific ex- 
perimental conditions turn out to be very critical for 
DNA transport measurements. Thus, experiments have 
yielded contradictory results as to the conduction prop- 
erties of DNA and are rather difficult to analyze. DNA 
has been characterized as a pure insulator ; 14 : 16 as a wide- 
band gap semiconductor ji 3 . and as a metallic system, 12 : 20 



Especially interesting are recent transport measurements 
on single poly(dG)-poly(dC) oligomers in aqueous so- 
lution, which displayed metallic-like I-V characteristics 
and an algebraic behavior in the length dependence of 
the conductance^ 

Notwithstanding this variety of results and the problems 
related to the experimental set-up, the possibility of us- 
ing DNA in molecular electronics is extremely attrac- 
tive since it would open a vast range of potential ap- 
plications because of its self-assembling and recognition 
properties^ Alternatively, DNA can also be used as a 
template in molecular electronic devices ^LSSiS 
From a theoretical point of view, the knowledge 
of the electronic structure of the base-pairs, the 
sugar/phosphate mantle and their mutual interactions 
is required in order to clarify the transport processes 
that may be effective in DNA. First principle approaches 
are the most suitable tools for this goal. However, the 
huge complexity of this molecule makes ab initio cal- 
culations still very demanding, so that only few inves- 
tigations have been performed, mainly in well-stacked 
periodic structures, 24 i 25 i 26 i 27 i 28 i 29 i 3Q i 31 i 32 i 33 i 34 To compli- 
cate this picture, environmental effects such as the pres- 
ence of water molecules and counterions which sta- 
bilize the molecular structure make ab initio calcula- 
tions even more challenging pSS 7 . Hence, Hamiltonian 
modela 36 ! 37 ! 38 ! 39 ! 40 ! 41 ! 42 ! 43 ! 44 : 45 that isolate single factors 
affecting electron transport are still playing a significant 
role and can help to shed more light onto the above issues 
as well as guide first principle investigations. 
Recently, Cuniberti et a/4i proposed a minimal model 
Hamiltonian to explain the semiconducting behavior pre- 
viously observed by Porath et a/4£ in suspended short 
(up to 30 base-pairs) poly(dG)-poly(dC) molecules. Re- 
markably enough, this experiment was performed on sin- 
gle molecules, in contrast to most transport experiments 
involving bundles of molecules. Molecular systems like 
Poly(dG)-poly(dC) (or Poly(dA)-poly(dT)) are very at- 
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tractive from a theoretical standpoint since, being pe- 
riodic, band-like transport as a result of 7r-orbitals hy- 
bridization may be more efficient than in its strongly dis- 
ordered counterparts, e.g. in A-DNA. The above model^i 
mimics the electronic structure of the complex poly(dG)- 
poly(dC)-backbone system by a tight-binding chain to 
which side chains are attached. Electrons can hop along 
the central chain but not along the side chains. As a 
result a gap in the electronic spectrum opens. The gap 
is obviously temperature independent and the transmis- 
sion near the Fermi level would show a strong exponen- 
tial dependence due to the absence of electronic states to 
support transport. 

An immediate issue that arises is how stable this elec- 
tronic structure, i.e. two electronic bands separated by a 
gap, is against the influence of several factors which are 
known to play an important role in controlling charge 
propagation in DNA molecules, such as static and dy- 
namic disorder 47 i 48 i 49 i 5Q i 51 i 52 i 53 and environment, 26 ! 27 ! 40 
In particular, the environment can act as a source of de- 
coherence for a propagating electron (or hole) ^ it can 
induce structural fluctuations that support or restrict 
charge motion^ or it can introduce additional electronic 
states within the fundamental gapii 9 ^ 7 . As it has been 
demonstrated experimentally, a modification of the hu- 
midity causes variations of orders of magnitude in the 
conductivity of DNAi^^ 5 . Moreover, the recent single- 
molecule experiments of Xu et alm& suggest that the en- 
vironment may strongly modified the low-bias transport 
properties of DNA oligomers. 

In this paper we elaborate on the role played by the en- 
vironment by addressing signatures of the bath in the 
electronic transmission spectrum of the DNA wire in dif- 
ferent coupling regimes: the mean-field approximation as 
well as weak-coupling and strong-coupling limits. Antic- 
ipating some of our results, we find that the semicon- 
ducting gap closes on the mean-field level as a result of 
thermal fluctuations. In the weak-coupling limit, how- 
ever, the gap opens with increasing temperature. In both 
cases the electronic gap is an "intrinsic" property of the 
system. On the contrary, a bath-induced pseudo-gap is 
formed in the strong coupling limit, i.e. an energy region 
with a low (but finite) density of electronic states. We 
have further found in this regime that the transmission 
at the Fermi level exponentially decreases with the wire 
length L, t(E-p) ~ c lL . The decay rate 7 is however 
rather small ~ 0.2 A -1 . This together with a noticeable 
dependence of 7 on the electron-bath coupling clearly 
indicates that incoherent pathways do appreciably con- 
tribute to charge transport in the strong coupling limit. 

In the next section we introduce the model Hamiltonian 
and derive the corresponding Green functions which are 
required to calculate the linear conductance. In section 
III different approximation schemes associated with dif- 
ferent coupling regimes to the bath are discussed. The 
influence of structural disorder on our results is also pre- 
sented. Finally, our summary follows in section IV. 
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FIG. 1: (Color line) Schematic drawing of the DNA molecular 
wire in contact with a dissipative environment. The central 
chain with N sites is connected to semiinfinite left (L) and 
right (R) electronic reservoirs. The bath only interacts with 
the side chain sites (c), which we call backbone sitesi— The 
Hamiltonian associated with this model is given by Eqs. 0, 
@, and © in the main text. 



II. HAMILTONIAN MODEL 

Along the lines of Refs. l4ll we represent the DNA molecu- 
lar wire containing ./V base pairs by the following nearest- 
neighbour tight-binding Hamiltonian (sec Fig.0: 



Wei = ^E^-^iE^v 

3 3 

j 

3 

= Tic + 7~tb + %c 



H.c. 



(1) 



Hereby Tic and H-b are the Hamiltonians of the central 
and side chains, respectively, and TLc-h is the coupling 
between them, tii and t± are hopping integrals along the 
central chain and between the backbone sites and the 
central chain, respectively. If not stated otherwise, the 
on-site energies will be later set equal to zero. The He 
Hamiltonian can be considered as effectively modeling 
one of the frontier orbitals of the poly(dG)-poly(dC) sys- 
tem, e.g. the highest-occupied molecular orbital, which 
is localized on the guanine bascsi^ 7 ^ 8 . The side chain in- 
duces then a perturbation of the 7r-stack leading to the 
opening of a temperature independent semiconducting 
gap in the electronic spectrum, the gap being propor- 
tional to the transversal hopping integral t± 4i Since this 
model shows electron-hole symmetry, two electronic man- 
ifolds containing N states each, are symmetrically situ- 
ated around the Fermi level, which is taken as the zero 
of energy. 

We focus here on the influence of the environment on the 
electronic structure and consequently on the transport 
properties of the model described by Ti. e \- As it has been 
demonstrated in the past years, correlated fluctuations 
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of hydrated counterfoils strongly influence electron(holc) 
motion along the double-helixpi2i2& Recent Raman and 
neutron scattering experiments on lysozymc have shown 
that the protein dynamics follows the solvent dynam- 
ics over a broad temperature range. Especially, confor- 
mational changes, low-energy vibrational excitations and 
the corresponding temperature dependences turned out 
to be very sensitive to the solvents dynamics^ We con- 
sider the vibrational degrees of freedom of countcrions 
and hydration shells in DNA as a dynamical bath able to 
act as a dissipative environment. In this model Hamilto- 
nian approach, we do not consider specific features of the 
environment but represent it by a phonon bath of M har- 
monic oscillators. We further make the assumption that 
the bath is only directly affecting the side chain whereas 
the central chain is well screened by the latter. Then, the 
extended Hamiltonian becomes: 



= H c \ + Hb + H c 



(2) 



where Hb and H c ~b are the phonon bath Hamiltonian 
and the backbone-bath interaction, respectively. B a is 
a bath phonon operator and A Q denotes the electron- 
phonon coupling. Note that we assume a local coupling 
of the bath modes to the electronic density at the side 
chain. Later on, the thermodynamic limit (M — > oo) in 
the bath degrees of freedom will be carried out and the 
corresponding bath spectral density introduced, so that 
at this stage we do not need to further specify the set of 
bath frequencies Q a and coupling constants A Q . 

Finally, we include the coupling of the molecular wire to 
semiinfinite left (L) and right (R) electrodes: 



H = Hw + E eko-dk^rikfT 

keL,R,er 

+ E (Vk,i4^i+H. c .) 

keL, cr 

+ E (V KN di a b N + K.c.) 

kG-R,<r 

= Hw + Wl/R + Hh-C + Hr-C 



(3) 



The Hamiltonian of Eq. © is the starting point of 
our investigation. Performing the Lang-Firso\iS uni- 
tary transformation H = c s Hc~ s with the generator 
S = E a A X a /n a )c] Cj (B a - fit) and & = -S, the lin- 
ear coupling to the bath can be eliminated. In the re- 
sulting effective Hamiltonian only the backbone part is 
modified since the central chain operators be as well as 
the leads' operators are invariant with respect to the 
above transformation. The new Hamiltonian reads: 



H = He + Hl/r + Hb + ^l/r-c 

+ ( e -A)^c] Cj -txx;[^ 



H.c. 



(4) 



X = exp 



\2 

a = v ^. 



Let's define two kinds of retarded thermal Green func- 
tions related to the central chain Gjt {t) and to the back- 
bones Pjg(t), respectively (h = 1): 

G jt {t) = -»e(t)^ [&,(<), 6j(0)] + \, (5) 

p jt {t) = -ie(t)([ Cj (t)Af(t),cJ(o)^(o) 



where is the Heaviside function and the average is 
taken w.r.t. H. With the above definitions and using 
the equation of motion technique (see Appendix ^) we 
arrive to an expression for the Fourier transform of the 
central chain Green function which reads, to lowest-order 
in tj_: 



G-\E) = GoHE)-tiP(E) 

G^(E) = £1-W c -£l(£)-£r(£). 



(0) 



In this equation Gq(E) is the Green function of a chain 
without backbones and connected to the left and right 
electrodes. The influence of the latter is comprised in the 
complex self-energy functions El/r(£ , )i^ The polaronic 
Green function P(-E) is explicitly given by: 

/>OG 

P tj (E) = -i5 tj / dte^+^e- 1 
Jo 

x "(l-/c)c-* (t) +/ c e -*<-*>] (7) 

withe-*W = (X(t)X^{0)) B being a dynamical bath cor- 
relation function. The average (-) B is performed over the 
bath degrees of freedom. Working to lowest order in t± 
allows to use a zero-order Green function for the side 
chain in Eq. (7), i.e. G c oij (t) ~ <%e ( e ~ A > *. f c is the 
Fermi function at the backbone sites. In what follows 
we consider the case of empty sites by setting f c = 0. 
Note that P is a diagonal matrix, i.e. it only modifies 
the on-site energies in the Hamiltonian. 
In order to get closed expressions for the bath thermal 
averages it is appropriate to introduce a bath spectral 
density 5 - defined by : 



j{uj) = V \ 2 j{uj - n a ) = j (-) s c -^eH, (8) 



where uj c is a cut-off frequency related to the bath mem- 



ory time r c 



It is easy to show that the limit uo c 



oo corresponds to a Markovian bath, i.e. J(t) ~ JaS(t). 
Using this Ansatz, $>(t) can be written in the usual way»£& 



J(u) 



1 



^ 1 — cos ivt 

' eP" - 1 



(9) 



Although the integral can be performed analytically^, 
we consider in some limiting cases where it is easier 
to work directly with Eq. |JS}- 

In the transport calculations, we limit ourselves to treat 
the low voltage regime, thus neglecting non-equilibrium 
effects as well as the inelastic part of the total current. 
As a result, one can still define a linear conductance g as 
follows i^i 



9(E) 



h J \ dE 
t(E) = Tr{A L GA R Gt} 



t(E), 



(10) 



where Al,r = i ^£l,r — £ L r] are spectral densities of 
the leads. Although the foregoing expression is similar 
to Landaucr's formula, we stress that the influence of 
the phonon bath does implicitly appear via the Green 
function G. Hence, both coherent and incoherent path- 
ways for charge transport mediated by phonon processes 
are included in Eq. (10). We concentrate our discussion 
on the temperature and length dependence of t(E). In 
what follows we always plot t(E) rather than g to filter 
out temperature effects arising from the derivative of the 
Fermi function in Eq. (10). For completeness the current 
as given by I(V) = (2e/h) J dE (f(E - eV/2) - f(E + 
eV /2))t(E) is also shown. We remark however, that this 
expression neglects non-equilibrium effects, which are be- 
yond the scope of this investigation. 



III. LIMITING CASES 

We use now the results of the foregoing section to dis- 
cuss the electronic transport properties of our model in 
some limiting cases for which analytic expressions can 
be derived. In all cases, we use the wide-band limit 
in the electrode selfenergies, i.e. £l,£?(-E) = — i^h^U^ij 
and Er^j(S) = — iTjiSNiSjsrj. We discuss the mean- 
field approximation and the weak-coupling regime in the 
electron-bath interaction as well as the strong-coupling 
limit. Farther, the cases of ohmic (s = 1) and super- 
ohmic (s = 3) spectral densities are treated. 



A. Mean-field approximation (MFA) 

Within the mean-field approximation bath fluctuations 
contained in P(E) are neglected. The MFA can be intro- 
duced by writing the phonon operator X as (X) B + SX in 
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FIG. 2: (Color online) Electronic transmission and corre- 
sponding current in the mean-field approximation for two 
different temperatures. Parameters: N = 20, Jo/W = 
0.12,tj./t|| =0.5,r L/R /t|| =0.5. 



0(SX). As a result a real, static and temperature depen- 
dent term in Eq. JHJ is found: 



G-\E) = G^(E)-tl 



E-e + A + i0+ 



where |(<-f) B | 



-2«(t) an( j K (rp^ ig given 



k(T) 



[ ^J(oj)coth 
Jo ^ 



2k B T 



(11) 



(12) 



The effect of the MF term is thus to scale the bare 
transversal hopping t± by the exponential temperature 
dependent factor e -K ( T ). 

In the case of an ohmic bath, s = 1, the integrand in 
k(T) scales as 1/lo p , p = 1, 2 and has thus a logarithmic 
divergence at the lower integration limit, see Eqs. JSJ and 
(12). Thus, the MF contribution would vanish. In other 
words, no gap would exist on this approximation level. 
In the superohmic case (s = 3) all integrals are regular. 
One obtains A = J du; tu^ 1 J(uj) = T(s — 1) J = 2J , with 
r(s) being the Gamma function and k(T) reads: 



2J„ 



2(^) 2 Ch(2, 



k B T 



- 1 



(13) 



(,n(s,z) = X)^Lo( n Z )~ S IS ^ ne Hrrrwitz ^-function, a 
generalization of the Ricmann £-function*£ 
It follows from Eq. (13) that k(T) behaves like a con- 
stant for low temperatures (feT/wc < 1), k(T) ~ Jo/w c , 
while it scales linear with T in the high-temperature limit 
(k B T/uj c > 1), k{T) ~ J /u c (l + 2k B T/u c ). 
If Jq vanishes, A is zero and (X) B = 1- Thus we re- 
cover the original model of Ref. |4l| which has a gap pro- 
portional to t±. For Jo ^ and at zero temperature 

_Jo 

the hopping integral is roughly reduced to t±e w = which 
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FIG. 3: (Color line) Electronic transmission and correspond- 
ing current in the weak-coupling limit with ohmic dissipa- 
tion (s = 1) in the bath. Parameters: N = 20, Jo/u> c = 
0.2, t± A,, =0.6,r L /R/t|| =0.5 



is similar to the renormalization of the hopping in Hol- 
stein's polaron model^ though here it is t± rather than 
fn the term that is rescaled. At high temperatures t± 
is further reduced (k(T) ~ T) so that the gap in the 
electronic spectrum finally collapses and the system be- 
comes metallic, see Fig. [21 An appreciable temperature 
dependence can only be observed in the limit Jq/oj c < 1; 
otherwise the gap would collapse already at zero tem- 
perature due to the exponential dependence on J . We 
further remark that the MFA may be only valid in the 
regime Jq/lu c < I, kp,T/uj c < 1, otherwise multiphonon 
processes in the bath, which are not considered at this 
stage, become increasingly relevant. 



B. Beyond MF: weak-coupling limit 

As a hrst step beyond the mean-field approach let's first 
consider the weak-coupling limit in P(E). For Jq/oj c < 1 
and not too high temperatures (k&T/u) c < I) the main 
contribution to the integral in Eq. J2J) comes from long 
times t ^> uj^ 1 . With the change of variables z = cut, 



<!>(£) can be written as: 



S(t) - J LO~ S t 




x 1 -e 



As far as u c t ^> f3uj c this can be simplified to: 



*(t) R» J Q uj- s t L - s I dxr 3 - 2 ' 



(14) 



x 1 - e 



( puj c 1 — cos z 

UJ r t z 



(15) 



Since in the long-time limit the low-frequency bath modes 
are giving the most important contribution we may ex- 
pect some qualitative differences in the ohmic and super- 



ohmic regimes. For s = 1 wc obtain $(t) 
which leads to (using A(s = 1) = Jo): 



. Jo k-&T 



G-\E) = G-\E)-tl 



1 



E+J +iir^k B T 



;1.(16) 



i.e. there is only a pure imaginary contribution from the 
bath. For the simple case of a single site coupled to a 
backbone one can easily see that the gap approximately 
scales as V^bT; thus it grows with increasing temper- 
ature. This is shown in Fig. [31 where we also see that 
the intensity of the transmission resonances strongly goes 
down with increasing temperature. The gap enhance- 
ment is induced by the suppression of the transmission 
peaks of the frontier orbitals, i. c. those closest to the 
Fermi energy. 

For s = 3 and k^TjLo c < 1, $(t) takes a nearly tem- 
perature independent value proportional to Jq/uj c . As 



a result the gap is slightly reduced (tj 



1 1 e 



but, because of the weak-coupling condition, the effect is 
rather small. 

From this discussion we can conclude that in the weak- 
coupling limit ohmic dissipation in the bath induces an 
enhancement of the electronic gap while supcrohmic dis- 
sipation does not appreciably affect it. In the high- 
temperature limit k^T/ujc > 1 a short-time expansion 
can be performed which yields similar results to those 
of the strong-coupling limit (see next section)^ so that 
we do not need to discuss them here. Note farther that 
the gap obtained in the weak-coupling and mean-field 
limits is an "intrinsic" property of the electronic sys- 
tem; it is only quantitatively modified by the interac- 
tion with the bath degrees of freedom. We thus trivially 
expect a strong exponential dependence of t(E = E-p) 
on the wire length, typical of virtual tunneling through 
a gap. Indeed, we find t(E = Ep) ~ exp (— (3 L) with 



C. Beyond MF: strong coupling limit (SCL) 



In this section we discuss the strong-coupling regime, as 
defined by the condition Jo/lu c > 1. This may be the 
regime to be found in presence of an aqueous environ- 
ment, as recent theoretical estimations using the classi- 
cal Onsager model for solvation processes have shown^ 
In the SCL the main contribution to the time integral in 
Eq. arises from short times. Hence a short-time ex- 
pansion of <&(t) may already give reasonable results and 
it allows, additionally, to find an analytical expression for 
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FIG. 4: (Color line) Temperature dependence of the real and 
imaginary parts of P(E) for N = 20, Jo/W = 10,tj_/i|| = 
0.4, r L / R /t|| = 0.5. With increasing temperature the slope of 
the real part near E = decreases and the imaginary part 
broadens and loses intensity. A similar qualitative dependence 
on Jo was found (not shown). 



P{E). At t < cj" 1 we find, 



*(*) 



iAt+(u c t) 2 k (T) 
-iS ej 



(17) 



= -i* 







dte i(B_e+i0+) * e" ( " c * )2Ko(T) 

(E-e + iO+) 2 



1 



exp 



1 + erf 



u c \Jku(T) 

i(£-e + i0+) 



2aj c y / K (T) 



«o(T) 



dw JUjS) coth 



Before presenting the results for the electronic transmis- 
sion, it is useful to first consider the dependence of the 
real and imaginary parts of P{E) on temperature and on 
the reduced coupling constant Jo/u> c - Both functions are 
shown in Fig. 0] We see that around the Fermi level at 
E = the real part is approximately linear, Re P(E) ~ E 
while the imaginary part shows a Gaussian-like behavior. 
The imaginary part loses intensity and becomes broad- 
ened with increasing temperature or Jo, while the slope 
in the real part decreases when /cb? 1 or Jo are increased. 
If we neglect for the moment the imaginary part (the 
dissipativc influence of the bath) , we can understand the 
consequences of the real part being nonzero around the 
Fermi energy, i.e. in the gap region. The solutions of the 
non-linear equation det\(E — t 2 ± Re P(E))1 — Hc\ = give 
the new poles of the Green function of the system in pres- 
ence of the phonon bath. For comparison, the equation 
determining the eigenstates without the bath is simply 
det\(E - t\/E)l - Hc\ = 0. It is just the l/E depen- 
dence near E = that induces the appearance of two 
electronic bands of states separated by a gap4i In our 
present study, however, Re P(E — > 0) has no singular be- 
havior and additional poles of the Green function may be 




FIG. 5: (Color line) Upper panel: t(E) with ImP(E) = 0; 
the intensity of the resonances on the central narrow band is 
strongly dependent on Jo/u c and k&T (not shown). Temper- 
ature dependence of t(E) with full inclusion of P(E) (middle 
panel) and corresponding current (lower panel) for N = 20, 
Jo/Wc = 5,tx/tii = 0.5,r L / R /t|| = 0.2. The pseudo-gap in- 
creases with temperature. 



expected to appear in the low-energy sector. This is in- 
deed the case, as shown in Fig. [5] (upper panel). We find 
a third band of states around the Fermi energy, which we 
call a polaronic band because it results from the strong 
interaction between an electron and the bath modes. The 
intensity of this band as well as its band width strongly 
depend on temperature and on Jo. When k^T (or Jo) 
become large enough, these states spread out and even- 
tually merge with the two other side bands. This would 
result in a transmission spectrum similar of a metallic 
system. 

This picture is nevertheless not complete since the imag- 
inary component of P(E) has been neglected. Its inclu- 
sion leads to a dramatic modification of the spectrum, 
as shown in Fig. [S] (middle panel). We now only see 
two bands separated by a gap which basically resembles 
the semiconducting-type behavior of the original model. 
The origin of this gap or rather pseudo-gap (see below) 
is however quite different. It turns out that the imagi- 
nary part of P{E), being peaked around E = 0, strongly 
suppresses the transmission resonances belonging to the 
central band. Additionally, the frontier orbitals on the 
side bands, i.e. orbitals closest to the gap region, are 
also strongly damped, this effect becoming stronger with 
increasing temperature (ImP(£) broadens). This latter 
effect has some similarities with the previously discussed 
weak-coupling regime. Note, however, that the new elec- 
tronic manifold around the Fermi energy does not ap- 
pear in the weak-coupling regime. We further stress that 
the density of states around the Fermi level is not ex- 
actly zero (hence the term pseudo-gap); the states on 
the polaronic manifold, although strongly damped, con- 
tribute nevertheless with a finite, temperature dependent 
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incoherent background to the transmission. As a result, 
with increasing temperature, a crossover from tunneling 
to activated behavior in the low- voltage region of the I-V 
characteristics takes place, see Fig. [5] (lower panel). The 
slope in the I-V plot becomes larger when t± is reduced, 
since the side bands approach each other and the effect 
of ImP(E) is reinforced. 

In Fig. an Arrhenius plot of the transmission at the 
Fermi energy is shown for different strengths of the 
transversal hopping integral and the electron-bath cou- 
pling. After a nearly T-independent region, the transmis- 
sion strongly grows up following approximately a e ~ 1 ' T 
law. Increasing the coupling to the phonon bath makes 
the suppression of the polaronic band around E = less 
effective (lmP(E ~ 0) decreases) so that the density 
of states around this energy becomes larger. Hence the 
absolute value of the transmission also increases. Simi- 
lar T-dcpcndcnccs have been experimentally observed in 
poly(dG)-poly(dC)i£ as well as in A-DNA.il On the other 
side, increasing ij_ leads to a reduction of the transmis- 
sion at the Fermi level, since the energetic separation of 
the side bands increases with tj_. 

We have further investigated the length dependence of 
the transmission at the Fermi energy. This is a very 
important aspect that helps to identify the influence of 
different transport mechanisms HJik The results are dis- 
played in Fig. [7| for different values of the reduced cou- 
pling Jq/lu c . For a homogeneous chain (on-site ener- 
gies are set to zero) an exponential dependence on the 
chain length t(Ep) ~ e~~ tL was found. In this expres- 
sion L = 7Vao, where N is the number of sites on the 
molecular wire and ao ~ 3.4 A is the average distance 
between consecutive base pairs. Note that the inverse 
decay lengths 7 are rather small ~ 0.1 — 0.3 A . An 
exponential dependence usually indicates virtual tunnel- 
ing through a gap. Inverse decay lengths, as extracted 
c. g. from complex band structure calculations 



however much larger that those obtained in the present 
investigation. So have recent DFT-bascd calculations 
found values of 7 ~ 1.5 A for gap tunneling in dry 
Poly(dG)-Poly(dC) oligomers^ With increasing bath 
coupling the exponential dependence farther weakens and 
eventually becomes algebraic t(Ep) ~ N~ a . The intro- 
duction of a tunnel barrier as realized e.g. through in- 
sertion of (AT)„ groups, by shifting the on-site energies 
along a finite segment of the chain increases the inverse 
decay length 7 by a factor of 2, approximately. Obvi- 
ously, this model cannot describe the crossover from su- 
perexchange mediated electron transfer (strong exponen- 
tial behavior) to sequential hopping-mediated transport 
(algebraic dependence) as a function of the wire length 
N, as discussed in other worksiii*^ We guess that vibra- 
tional excitations inside the central chain, which renor- 
malize the longitudinal hopping integral to, have to be 
included to get this non-monotonic transition. 
From the previous discussion we may conclude that elec- 
tron transport on the low-energy sector of the trans- 
mission spectrum is supported by the formation of po- 
laronic states. Though strongly damped, these states 
manifest nonetheless with a finite density of states inside 
the bandgap. 

It has been meanwhile demonstrated 26 i 48 i 49 ' 50 i 51 i 52 i 53 
that electron (or hole) motion in DNA is extremely sensi- 
tive to different kinds of disorder: static disorder (random 
base-pair sequences), structural fluctuations and inhomo- 
geneities of the counterions distribution along the back- 
bones. These factors may strongly distort the base pair 
stacking along the double helix and eventually affect the 
electronic transport properties. They deserve a separate 
study. However, as a test for the stability of our results 
we have randomly varied the on-site energies along the 
central chain by extracting them from a Gaussian distri- 
bution with variance ctq. In this way we are simulating 
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that a bath-induced pseudo-gap in the electronic spec- 
trum can appear for strong electron-bath coupling giving 
a temperature-dependent background around the Fermi 
energy. As a result the system may show with increasing 
temperature a transition from a tunneling to an activated 
behavior in the low-bias region when coupled to an ex- 
ternal dissipative bath. An Arrhenius-likc temperature 
dependence of the transmission at the Fermi level and a 
rather weak exponential dependence on the wire length 
were additionally found, indicating a strong contribution 
of incoherent pathways of the charge carriers. 
A natural extension of this investigation would be the 
inclusion of non-equilibrium effects at large bias and 
consequently of inelastic components of the current. 
This issue is although interesting from a formal point of 
view, since the Lang-Firsov transformation introduces 
polaronic rather than pure electronic propagators, 
see Eq. (5). For the former the appropriate Keldysh 
Green functions should be derived in order to deal with 
the non-equilibrium regime. This problem deserves a 
separate investigation which is now in progress. 



fluctuations inside the central chain. In Fig. 8 the cases 
of weak (cto ~ 0.12tn) and strong disorder (<jq ~ to) are 
shown. Two main features can be seen: (i) the transmis- 
sion resonances on the side bands are strongly washed out 
and lose in intensity, and (ii) the pseudo-gap is slightly 
reduced with increasing disorder. However, the suppres- 
sion of the central band due to ImP(E) and hence, the 
pseudo-gap formation is not affected by this kind of dis- 
order. As soon as electronic states shift from the side 
bands into the region with nonzero Im P(E) they are 
strongly damped and thus the pseudo-gap structure of 
the spectrum is conserved. A similar effect of disorder 
is expected in the other coupling regimes to the bath 
degrees of freedom discussed above. 



IV. SUMMARY 

Charge propagation in DNA molecules is extremely sen- 
sitive to disorder and environmental effects. We have 
focused in this paper on the influence of a dissipative 
environment on the electronic transport properties of a 
model Hamiltonian which mimics some basic features of 
the electronic structure of DNA oligomers. Although we 
have chosen Poly(dG)-Poly(dC) molecules as a reference 
point, we believe that our model is quite generic and may 
be useful for a large class of 7r-conjugated systems. 
We have shown that a mean-field approximation can- 
not fully catch the action of a dissipative environment 
on charge transport, because it only gives a real, energy 
independent contribution. Indeed, while the mean-field 
approach leads to gap reduction with increasing tem- 
perature, bath fluctuations eventually lead to gap open- 
ing in the weak-coupling limit. We have further shown 
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APPENDIX A: DERIVATION OF EQ. © 

The equation of motion for the retarded Green function 
in Eq. (JSJ in the frequency representation reads: 

EG ej {E) = ([b j ,b i } + ) + {{[b j ,n] 

Using it we get for the Hamiltonian of Eq. (J3J : 

£ [Go\E)] en G n j (E) = <% - t ± {{ctX\bt))(Al) 

n 

[Gq 1 ^)! = (E-e b )5 n t+t\\(Sn,l+l+SnJ-l) 
L J in 

— £lAl<5„l — T,B,5iNS n N 
„ V"^ |Mc,l(JV)| 2 

~ 2-f £- ek+ i0+ 

Now, equations of motion from the "right" may be writ- 
ten for the Green function Zf } (E) = ((c^Afjftt)), 
leading to : 

J2zL(E)[G^(E)] m . = - tx (te*| C t*t)) 

rn 

= -t x P tj (E) (A2) 
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In the former equations we have neglected cross-terms of 
the form ((ceX\cj)), since they will give contribution of 
0(t 3 j_). Inserting Eq. (|A2I> into Eq. (|A1|) we arrive at the 
matrix equation: 

G(E) = G (E) + G (E)Z B (E)G (E), 

which can be transformed into a Dyson-like equa- 
tion when introducing the irreducible part 'Sb(E) = 



^W(E) + ^(E)G (E)^(E) + ...: 

G(E) = G (E) + G (E)^(E)G(E). (A3) 

From Eq. (|A3|) it immediately follows Eq. © with 
Sg r (i?) = t 2 ± P(E). We emphasize that these expres- 
sions are exact only to lowest-order in the transversal 
hopping t^. This approximation may be justified in the 
low-voltage limit we are dealing with. 
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